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Abstract 

We construct a model for the universe based on the existence of quantum fields at finite tem- 
perature in the background of Robertson- Walker spacetime in presence of a non-zero cosmological 
constant. We discuss the vacuum regime in the light of the results obtained through previous 
studies of the back-reaction of massless quantum fields in the static Einstein universe, and we 
argue that an adiabatic vacuum state and thermal equilibrium is achieved throughout this regime. 
Results shows that such a model can explain many features of the early universe as well as the 
present universe. The model is free from the basic problems of the standard Friedmann cosmology, 
and is non-singular but involves a continuous creation of energy at a rate proportional to the size 
of the universe, which is lower than that suggested by the steady-state cosmology. 
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I. INTRODUCTION 

The Friedmann models for the universe which were deduced from the Einstein field 
equations with a vanishing cosmological constant and vanishing covariant derivative of the 
energy-momentum tensor, described an expanding universe that starts from a singularity at 
t = O.This develops into a universe that is either: (a) ever expanding with acceleration if the 
curvature is negative (k = —1), or (b) ever expanding with an ultimate constant speed if the 
geometry of the spatial section is flat (k = 0), or (c) expanding until reaching a maximum 
radius, then collapsing under its own gravitational field to an ultimate singularity, a model 
which is described as a closed spacetime with positive curvature (k = 1). These models 
found their applications in two main observational achievements of cosmology during of the 
last century; the Hubble discovery of the redshift-distance relationship, and the discovery of 
the Cosmic Microwave Background radiation (CMB). 

The existence of a homogenous and isotropic CMB was a prediction of the Gamow and col- 
laborators's big bang theory, which has assumed that the universe started as a homogeneous 
and isotropic distribution of particles in thermal equilibrium at very high temperatures. As 
the universe expanded and cooled down, massive particles decoupled from thermal equilib- 



2 



rium. Accordingly radiation was released as electrons combined with hydrogen and helium 
nuclei to form atoms. These radiations are believed to have left a relic which is the detected 
CMB. This scenario was further developed and refined and was called the Standard Big 
Bang model (SBB). 

The SBB utilized statistical physics, particle physics and the standard general relativity 
in which particles were not considered in their field theoretic description, but rather were 
dealt with through their general phenomenological description. On the other hand the 
SBB model did not start with the universe at t = since physics cease to work in such a 
limit. For these reasons some basic questions remained unanswered; the existence of the 
initial singularity posed a challenge for the physicists to define the initial conditions for the 
universe. Indeed it was only assumed that the universe originated from vacuum without 
giving any explanation for such a birth. Quantum effects were not investigated within the 
standard model and therefore, some important contribution is missing from the SBB model. 

Recent analysis of the main features of the CMB suggested that the universe is highly 
homogeneous and isotropic on very large scale and that, if it is to be described geometrically 
by the Friedmann solutions, then it has to be nearly flat (i.e. k = 0) Lp. This result was 
confirmed by subsequent observations from WMAP by Spergel et al. [2[ and also later by 
Dunkley et al. [3J]. This conclusion was based on the fact that direct observational data 
suggests that the total matter and radiation density in the universe is about the same as 
the critical density needed to flatten the universe. On the other hand, it well known that 
the SBB model suffers from some basic shortcomings, these are: the horizon, the flatness, 
the magnetic monopoles and the formation of large cosmic structures problems. To resolve 
these problems Guth [4( was the first to suggest that the universe may have experienced 
a state of inflation during very early times of its development, such that the flat geometry 
of the universe was attained at a very early stage. This theory was further developed into 
a main trend in cosmologv where a diversity of inflationary models were suggested (for a 
comprehensive review see [5|). A model of inflation typically amounts to choosing a form for 
the potential, perhaps supplemented with a mechanism for bringing inflation to an end, and 
perhaps may involve more than one scalar field. In an ideal world the potential would be 
predicted from fundamental particle physics, but unfortunately there are many proposals for 
possible forms. It has become customary to assume that the potential can be freely chosen, 
and then one seek to constrain it with observations. 
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However, some authors 6| argue that explanations provided by inflation for the homo- 
geneity, isotropy, and flatness of our universe are not satisfactory, and that a proper expla- 
nation of these features will require a much deeper understanding of the initial state of our 
universe. On the other hand, and although inflationary models are spectacularly successful 
in providing an explanation of the deviations from homogeneity, these authors point out that 
the fundamental mechanism responsible for providing deviations from homogeneity, namely, 
the evolutionary behavior of quantum modes with wavelength larger than the Hubble radius, 
will operate whether or not inflation itself occurs. However, if inflation did not occur, one 
must directly confront the issue of the initial state of modes whose wavelength was larger 
than the Hubble radius at the time at which they were born. Under some simple hypotheses 
concerning the "birth time" and initial state of these modes (but without any fine tuning), it 
is shown that non-inflationary fluid models, in the extremely early universe, would result in 
the same density perturbation spectrum and amplitude as in inflationary models, although 
there would be no "slow roll" enhancement of the scalar modes f|. Other authors believe 
that inflationary theories are incomplete since it does not deal with basic puzzles such as 



the initial singularity, nor with dark energy indicated by the recent observations [7]. 

Ozer and Taha 8| devised a model universe free from the basic problems of the SBB model 
by assuming a universe predominantly kept at the state of critical density and, accordingly 
the scale factor of their model is given by 



R 2 (t) = R 2 + t 2 . (1) 

On the other hand Chen-Wu |9( adopted the prescription that the cosmological constant 
varies like l/a 2 (t). Their study resulted in a model that has the same form for the variation 
of the scale factor as in the SBB model, therefore would have no problem integrating the 
Gamow explanation of natural abundance, but it would again lead to a singular universe this 
time with continuous particle creation at a rate comparable to that suggested by the steady 
state theory. The Chen-Wu suggestion was shown to explain the cosmological constant 
problem through a phenomenological approach. The good feature of the Ozer and Taha 
model is the fluent removal of the SBB problems without the need for an inflation stage. 
However it remains that the assumption of having a universe starting up with a density 
exactly equal to the critical density will surely need justification. 

It is quite possible that the scale factor do not follow a monotonic behavior during the 



whole history of the universe, rather it is quite expectable that it follows different schemes 
during the different stages of the development of the universe. It is certain now that the 
universe has passed through many phase-transition states, that had different variations of 
the scale factor. If so, then one can say that the Friedmann models will stand as a simplified 
picture for the development of the universe and that it is good for describing the basic theme 
only. The generally accepted scenario now stipulates that there are at least three regimes: 
the vacuum-dominated regime which can be described by the control of the cosmological 
constant, the radiation-dominated regime during which massive particles were in an ultra- 
relativistic state in thermal equilibrium with radiation, and the matter-dominated regime 
which developed when particles attained non-relativistic states and settled to form a dust- 
like fluid. 

In this paper we construct a model for the universe with the background geometry of 
Robert son- Walker metric endowed with quantum fields at finite temperature. We assume 
the presence of a non-zero cosmological constant and seek the solution of the Einstein field 
equation. Therefore, the present model does not satisfy the standard Friedmann paradigm, 
and in order to construct the model we utilize the results of previous studies and calculations 
of the vacuum energy density in the static Einstein universe. The universe is shown to have 
a violent start from a non-singular Plank sized patch developing through the interaction 
of vacuum energy and curvature into a very hot spot and then transiting smoothly into 
a thermal universe that coasts into the present one. The model is free from the standard 
problems of the SBB based on the standard closed Friedmann model including the initial 
singularity which get smoothed-out by the quantum- vacuum effects. 



II. THE STATIC EINSTEIN UNIVERSE 



The investigation of the back-reaction effect of massless quantum fields at finite tempera- 
tures in the background of the static Einstein universe resulted in a relationship between the 
radius of the Einstein universe and its temperature showing some interesting features 11]. 
First, all solutions were shown to posess two regimes, the vacuum (Casimir) regime, through 
which the temperature rises from zero sharply reaching a maximum value of order of 10 32 i^ 
within a very small change of radii, and a Planck regime, through which the temperature 
decays exponentially to zero following a Planckian behavior (see Fig. 1). Second, it was 



5 



shown that the universe has no singular state; no static Einstein universe can be singular, 
rather the ground state is at zero temperature with a non-zero radius of order of Planck 
length. Third, a background (Tolman) temperature was calculated for states of Einstein 
universe with large radii, and it was found that the background temperature equals the 
observed value of 2.73 K at a radius nearly two order of magnitudes as compared to the 
present Hubble length. 

It is well known the static Einstein universe is a solution of the field equations with a 
cosmological constant, that can always be adjusted to balance the gravitational attraction 
of matter contained within the spatial section of the universe. The study of the variation of 
the value of the cosmological constant with temperature for successive states of the Einstein 
universe resulted in showing that the cosmological constant has large and nearly constant 
value during the Casimir regime, decaying according to the inverse square law (A ~ 1/a 2 ) 
during the Planck regime (see Fig. 2). The exceptional case was with the minimally coupled 
massless scalar field, where the cosmological constant was shown to have an infinite value 
at start, decaying exponentially to small values already within the Casimir regime 121 ] . 

III. A TIME-DEPENDENT COSMOLOGICAL CONSTANT 

The above mentioned results are quite motivating to consider a more realistic case of a 
time-dependent model in which the cosmological constant is changing with the scale factor 
according to the inverse square law. 

In order to include any and all sources of energy that would contribute to the total energy 
of the universe, let us write the Friedmann equations in the following form 



where (p m .r)eff is the energy density for mater and radiation, and Xcs represents any and 
all contributions coming from the cosmological constant or any other source of energy- 
momentum density including the vacuum energy density. In other form the Friedmann 
equation can be written as 




(2) 




(3) 
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where (To)totai = P m > r + Pv,with p v = (0|Tg \0) tot is the total vacuum energy density. Eq. 
can be re-written as 
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Now if 
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Comparing the left hand side of (J7J) with (J2]) we get 



% = 32tt<0|T |0> 



tot 



(5) 



(6) 



(7) 



(8) 



or k = 1 this will give the same result as that obtained in the case of closed static universe 



[Hi- 



From the basic Einstein field equation and for the Robertson- Walker metric we can deduce 
that 



T 0;, = P + 3 {^) (P + P)- 



(9) 



Applying the covariant derivative to the Einstein field equation we obtain 



/; : >(-) (P + P) = -^- 

a J 8n 



(10) 



This equation will stand as a replacement for the equation of state used by the standard 
Friedmann models where the right-hand side is taken to vanish on the assumption that A is 
constant. 

The second law of thermodynamics requires that 
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dE + pdV = TdS, 

which means that 



dE dV _ T dS_ 
dt ^ dt dt' 



(11) 



where E = m c 2 = (pV) , and V = 2ir 2 a 3 is the volume of the closed universe. Therefore 



dE 



— = Q n 2 a 2 p + 2ttV-^ 



da 



2^3 



dp 



dt 



dt 



so that for Eq. (TTTT) we obtain 



V 



p + 3 - )(p + p) 



dt 



T dS 
dt' 



(12) 



(13) 



where the dot denotes differentiation with respect to time. Substituting (Tl3]) in ffTO]) we get 



V 

TdS = dX. 

8vr 



(14) 



Therefore, a variable A may solve the entropy problem without the introduction of specific 
fields or irreversible processes. The idea that A may be variable has previously been suggested 



by several authors 



12| and [13]. 



For the homogeneous isotropic model with Robertson- Walker metric with p = p(t) 
p = p (t) and A = A (t) , Eq. ©, © and ® yield 



(15) 



']) ="" 1 (Pm, r (t)+A(t))-^, 

where a = (3/87r) and A = (A/871") . 

Thus dS/dt > requires that dA/dt < for all t > 0. Then From Eq.l flBll we may 
interpret A (t) as a vacuum or cosmological energy density. 



Now from ffTTl) and (fl4l we get 



djpV) dV V_ dX _ Q 
dt ^ dt 87r dt 



which means that 
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dt dt dt 

Taking p = \p for radiation filled universe we get 



d(pa 3 ) da 3 ndh , , 

+ P-3T + « 3 -37 = 0- ("I 



^ + ^ + ^ = 0. (17) 
da da a 



To solve this equation we consider the general form 

dy 

where the solution is given by 



(Jx + P(x)y = Q(x), (18) 



y = exp [— I (x)] J Q (x) exp [I (x)] dx + frexp [— I (x)] , (19) 



where 



I(x) = J P{x)dx, (20) 
and b is a constant. Here we have P (a) = - so that 

I ( a ) = f 1 P (a) da == In ( — \ (21) 

Jao \ a 0/ 



Also 



„ . . dA. . , 

Q(a) = --. (22) 

consequently we obtain the solution 

'-*(?) 4 -® 4 £m*5*' (23) 

We shall take po and ao to be the values of the energy density and scale factor at t — 0. 

If the total energy of the universe, including the vacuum energy, is to be constant then 
we have to take dE = 0, which implies that 

TdS + a 3 dA = 0. (24) 

One observes that the conditions A < and a > imply that the integral in (|23j) is 
negative so that the decrease of A as a increases generates a positive contribution to p. Thus 



the cosmological energy density A is continually depleted and transformed into radiation 
energy density in accordance with fl23|) . This consolidates the conjecture of Altaie and 
Setare 1121 1 . 

n 

Following |8j we take the initial time t = to be the moment when a = 0. Then if p = 0, 
Eq. (1231) is still valid and the empty curved space-time metric which is governed by 

K ~ \tfR = -Stt [2?] (vac) , f = 0, (25) 

soon will generate a non-empty universe, i.e., one with non- vanishing radiation density p. 
One might thus take Eq. (1231) to imply that radiation is being created out of the space- 
time curvature all times. The physical picture one has is therefore the one in which the 
universe is continuously created by the unfolding of space-time curvature. This is different 



from the continuous creation in the Bondi-Gold-Hoyle [18(] model which was suggested in 
the context of steady-state cosmology. The steady-state model is also based on classical 
general relativity with modified T^ v . The modification is by the addition of a covariant term 
which was chosen so that p, p and H remain constant throughout the de Sitter expansion 
of the universe. Continuous creation and the absence of an initial singularity are features of 
steady-state cosmology that are shared by the present model. 

In the present formulation the cosmological energy density A is related to entropy by 

TdS + a 3 dA = 0. (26) 

This equation may in fact be interpreted as an expression of the constancy of the total 
entropy of the cosmos, substance and spacetime, i.e., 

dS + dS c = 0, (27) 

where 

a 3 

dS c = -dA, (28) 

is the change in the entropy of the curved spacetime. On the long run (t — > oo) this change 
flattens the spacetime. Under these conditions one would intuitively expect a decrease in 
the entropy S c since the number of degrees of freedom one might associate with a state of 
high curvature should be larger than those associated with an almost flat space. 
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We note that in a universe of pure radiation the choice of A completely determines the 
model, since A determines p by Eq. (1231) . leaving Eq. ffl5|) to be solved for a (t) . However 
when, in addition to radiation matter is present, extra assumptions are needed to uniquely 
determine the model. In this case p = p m + p r , where p m is the rest-mass energy density 
and p T is the energy density of radiation and relativistic matter (to be referred simply as 
radiation in what follows). Then Eq. (Tl6l) reads 

d (p m « 3 ) + d (p r a 3 ) + pda 3 = -a?dk. (29) 

A plausible assumption that may readily be made is that the processes responsible for 
changes in rest-mass energy density are, except for the matter creation period, much slower 
than those responsible for creation of radiation, (see jsj]); i.e., 

(30) 

We also make the assumption that the pressure of the universe, under these conditions, 
is caused by its radiation 

P = Pr= ^Pr- (31) 

Eq. (fTEj) then yield 

d (p r a 3 ) + -p r da 3 = -a 3 dA, (32) 
3 

and 



dt 



(p m a 3 ) 



dt 



(p r a 3 ) 



d (p m a 3 ) « 0, (33) 

implying that Eq. ([16]) may be assumed to be approximately valid for the radiation com- 
ponent in both the radiation-dominated as well as the matter-dominated eras. This may 
alternatively be interpreted in the sense that a substance is created, by the unfolding of 
curved spacetime, as massless radiation. Yet another interpretation could be that the change 
in entropy of non-relativistic matter is much less than that of radiation, i.e., AS m A5 r , 
since — a 3 dA is a measure of the total change in the entropy. 
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IV. THE MODEL 



We now deduce a particular function for A, thereby defining a specific cosmological model 
in the classical class. Our starting point is the observation that the present value of the 
energy density of the universe is close to its critical value p c . Theoretically the investigations 
of the back-reaction effects of quantum fields at finite temperatures indicates that the density 
of the universe was always fixed at the critical value 

. , 3H 2 3 faV 
Pc(t ) = -— = — - . 34 

87T 87T \ a I 

From Eq. (|T5l) we notice that the condition p = p c requires that 




lr < 35 > 



So that A(t) comes to be 



AM = £ (36) 

The conditions A < and a > then immediately give k > so that k = 1. This is 
a significant deviation from the standard model where p = p c implies k = 0. This is an 
important feature of the present model, therefore, we can confidently conclude that 

A (.) = £. (37) 

The model now is completely specified in its dynamical structure. The rest of the work 
will mostly run in similar fashion to that of Ozer and Taha with a difference by a factor of 2 in 
some equations. However, as for the initial start of the universe we here do not need to follow 
the assumption of Ozer and Taha but would rather resort to take the results of previous 
works which defined for us the initial radius of the universe as a result of self-consistency 
condition applied on the Einstein field equations for a given quantum field source. The 
provisions of an equation of state determines the physical content. In the following sections 
we present the results of the calculations with some outlines. 
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A. The very early universe 



In our model the initial moment t = has been chosen to coincide with the state of the 
universe when its energy content is specified by the presence of Casimir energy resulting 
from the high initial curvature. Self-consistency of the Einstein field equation requires the 
universe to have a non-zero radius at T = 0. This was already calculated in a previous 
works (see [l^]) There are, initially, no such excitations since a (0) = and the total energy 
is locked up in potential form in spacetime curvature. With p = 0, Eqs. ( |23|) and (1371) yield 

>-M x -$y (38) 

Note that the condition po = requires ao 7^ 0. This implies that all functions a, p, S, 
T are finite at t — (as well as for all finite t as will soon become clear) so that the initial 
singularity of the standard model does not exist. 

Eq. (I38I) with p = p c = a(~ ) 2 , may now be solved for a (t) giving 

ada ( 1\ , . . 

This can be easily integrated to give 

1/2 

(40) 

It is clear that the universe starts accelerating during the Casimir regime and then soon 
get to an ultimate speed when it reaches the velocity of light. In this model we have no 
inflation but a direct parametric dependence of the radius on time. Note that in this solution 
of the field equations a —>■ 00 as t — > 00 i.e., the model is continuously expanding although 
k = 1. This is due to the variable cosmological energy density A in Eq. 0351) which renders 
the characterization of the asymptotic behavior being A— dependent. Different choices of A, 
for fixed k, can give different types of asymptotic behavior. The intrinsic geometry, for fixed 
t, is governed by the parameter k and is independent of A. 

The time-dependence of all functions in the model is completely determined by Eqs. ( |38l) 
and OH)]) . For the radiation energy density we have 



a(t) 



t 2 

2 
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The radiation temperature T is assumed to be related to p by 



P(T) 



1 



tc 2 N (T) T 4 



(42) 
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The effective number of spin degree of freedom iV (T) at temperature T is given by 
N (T) = iVb (T)+|iVf (T) , where (T) and Nf (T) refer to bosons and fermions respectively. 
We take N (T) to be constant throughout the pure radiation era. From Eqs. ( T4T|) and (142]) 
we obtain 



This is qualitatively the same result that was obtained by Ozer and Taha. In terms of 
the radius a the dependence of T is given by 



Thus, according to this model the universe have a cold start since T = at t = 0. For small 
a Q , this need, however, not to be different from a hot universe since temperature increases 
rapidly within a time-scale of order ao- Fig. 3 shows the qualitative time development of 
the temperature of the universe T according to Eq. fT44|) . Also, from Eqs. (1371) and fT40|) we 
get Fig. 4, which illustrates a qualitative relationship between T and A. 

The maximum temperature T max is obtained at t = V2a and is given by 




(43) 




(44) 




(45) 



For t > a , Eqns. (gT]) and (g3]) gives 



a 



(46) 



P 



and 




(47) 
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These equations are to be compared to those of the standard model, namely 

a 



PSM 



(2t) 



2 ' 



(48) 



and 



SM 



30a 



1/4 



1 

2t 



nl/2 



(49) 



Thus for t > ao the values of the energy density and temperature attained at a time 
t in the standard model are attained at time 2t in the present model. For t ~ ao, Eqs. 
(j4ip and fj4*3|) coincide with Eqs. (1481) and (j4*9]) of the standard model. One may, therefore, 
conclude that although this model is clearly different from the standard model in several 
aspects, such as having cold initiation, non-adiabaticity and regularity at t = 0, it possesses 
for T > T max essentially the same thermal history as the standard model. This is somewhat 
surprising since the dependence of a on t, in the present model, is completely different from 
that of the standard model: asM ~ t 1 ! 2 . On the other hand we should, however, expect 
some substantial deviation in the time-dependence as compared to the standard model if 
the cosmological problems of the standard model are to be avoided. 

If we have to look at the variation of the cosmological constant with the temperature 
of the universe, then we will obtain the temperature dependence depicted in Fig. 4. The 
figure shows a dependence which is similar to that we obtained for the Einstein universe for 
conformaly coupled scalar field. 

In particular, the dependence of the cosmic scale factor a on t determines the causal 
structure of the model. The horizon distance du (t) at time t is the proper distance travelled 
by light emitted at t = 



d>H (t) = a (t) 



dt 



(50) 



For the universe around us to be causally-connected to us at cosmic time t it is necessary 
that da (t) > ^proper (t) where d pr0 p er (t) is the proper distance, at time t between our galaxy 
and another galaxy most distant from us assuming that our galaxy is at r = 0. : 



^proper \t) Q (^) 



dr 



VI - kr 2 



a{t) 



'"max > k 



sinh V^, k = -1 

Sin Tmax; k = 1 



(51) 
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For k = and k — — 1 the universe is spatially infinite so that r max = oo. This implies 
that for k = and fc = — 1 global causal connection, i.e. causal connection for the whole 
space, is never established at any finite time t. The region of the universe which is causally 
connected at time t is limited in coordinate space to < r < r# (t) , where 



The horizon (or causality) problem for k = and k = — 1 in the standard model may be 
formulated only for the currently observed universe, i.e. in non-global terms. 

For k = 1 the universe is spatially finite and r max = 1. It is then possible to determine 
the time t = t caus when the whole universe is causally connected. This is given by 



For the standard model with k = 1, one finds that by the end of the radiation-dominated 
era, t r = 10 12 s, only a small part of the whole space is causally connected. To see this 



Using a = At 1 / 2 in the radiation-dominated era of the standard model where A = 
(2TT 2 N/15a) 1/4 a p T P . For N = 100, a p = 10 10 years and T P = 2.7K , the right-hand side 
of (PD is 0.014 so that Eq. ((521) with k = 1 gives r H (t r ) = sin 0.014 w 0.014 < 1. Global 
causal connection, i.e. up to r = 1, is realized during the matter-dominated era as may be 
verified using floU]) . Thus most of the currently observable universe was not in causal con- 
tact at the end of the radiation-dominated era. One is then unable to explain the observed 
isotropy of the background black-body radiation. 

In the present model Eqs. (j4"0l) and fl53|) determine the time t caus when global causality 
is established. One finds 




(52) 




(53) 




(54) 




(55) 



Now let 



t 



aosinhw =^ dt = y/2a (cosh w) dw, 



V2 
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then Eq. (1551) becomes 



V2a (cosh w) dw n 

(56) 



o (<2q + (a sinh w) ) 



2\i/2 2' 



tcaus = V2a sinh ( — ) = 1.9a . (57) 



which then would yield 



.2^2, 

Note that for the integral in (I53I) to converge it is necessary to have ao 7^ 0. Eq. (J5j 
indicates that global causal connection in the present model has been established at a very 
early time. Thus the present model does not possess a horizon problem. 

We observes that the "cold era" in this model is restricted to an interval when the 
causally connected part of the universe covers a tiny fraction of the whole space. When 
global causality is established at t caus = 1.9a the maximum temperature is surpassed and 
the whole universe attains the temperature T caus .From Eq. P5|) . we have: 



15a \ 1/4 / IN 1/4 



Ati 2 N / 

also with the help of Eqs. ( H5l) and ( 1571) . we get: 



(58) 



15a \ 1/4 / IN 1/4 



^ = 0-6916 ( — ) (3) . (59) 



this shows that 



T caus = 0.978T max (60) 



B. Radiation and matter 



The very early pure radiation era soon gives way to a period of matter generation. 
Throughout this period, a\ < a < ai say, Eq. ( |29|) is valid. For a > ai Eqs. (1321) and 
f!55|) hold, so that the total rest-mass energy 

E m = p m a 3 , (61) 
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remains approximately constant, i.e., 



Pn 



(62) 



Note that although the proper volume is 2ir 2 a 3 we shall, in accordance with general 
convention, take it to be simply a 3 . This is of no consequence since the measurable quantity 
is the energy density. 

The solution to Eq. (|32|) for p r with A given by Eq. (1371) is 



p = p exp 



In 



+ exp 



-In 



a E p 

_ | m 

a * A 



exp 



In 



(63) 



Now, if p = po when a = ao, then p becomes 



y 1 

a 



2^2 



2£ p 

m 



dd 



+ 



m 



1 +LU 



2\i 
P 
2 



+ 



m 



(64) 



so that 



m 

a 3 



2^2 



(65) 



Accordingly 



a 

2^2 



1 + cu 



a > a2, 



(66) 



where u is a dimensionless constant and a p is the present value of the scale factor. Note 
that although Eq. (|32|) for p r is the same as Eq. (TIB"}) , its solution (|66j) is not the analytic 
continuation of p in Eq. fl38|) . The reason is that the system is subject to different equations 
of state in the two regions: p = \p for ao < a < d\ and p = | (p— E^/a 3 ) for a > a 2 . These 
two regions do not, therefore, belong to the same phase. The region ai < a < a 2 during 
which rest-mass is created may, therefore, be thought of as a region of phase transition. This 
will be further discussed in the next section. 

When uj is expressed in terms of present values of p r and a one obtains 
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(67) 



The total radiation energy E T = p r a 3 is then 



E r 



a auj a r 
-a + 



a > a 2 . 



2 2 a 

For matter dominance to occur there must exist a value of a, a = a eq say, such that 



(68) 



E T (a eq ) = E^, 



(69) 



where for all a> a 2 , 



E r (a) ^ E^ when a ^ a, 



cq • 



(70) 



This requires that for a?, < a < a eq and at least in the neighborhood of a eq , a must be small 



enough for E T to be decreasing at a = a eq so that 



\ i *eq 
UJ > — - 



The value of a eq is given by 



(71) 



pl,al ± 



m p 



(pP a 3 ) 2 — 4a (2u;aap) 



1/2 



'-cq 



2a 



then from Eq. ([6 



(72) 



V2(p?) 1/2 



1/2 



(73) 



P^f« 3/2 (l+-) 3/2 /v/8(p?) 3/2 



"eq 



2a 



[(pP ) 2 a 3 (1 + a;) 3 /8 (p P ) 3 - 8a^a (1 + w) /2p?] 



1/2 



2a 



1 pPa 1 /2( 1+a; 



,3/2 



2(p?) 



px3/2 



1 ± 



1 - 



4a; 




l/2\ 


(1+co) 2 




J 



(74) 
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The temperature T eq at a = a eq is given by 



E? = pVa 3 p = \ 2 N(T P )T*a 



4„3 
P 



1 

30 



vr 2 iV (T cq ) T'a 



eq^eq 



Remembering that E T (a eq ) = E^ = p^a^, therefore Eq. (1751) becomes 



and 



so that 



with 



TpOp 



pi J \a. 



Now from Eq. (174 



a 1 / 2 (1+u) 
V2(p?) 



1/2 



pxl/2 



l eq 



2p\ 



1/2 



pPa 1/2 (l+^) 3/2 



2V8(p! 



PN3/2 



1 ± 



(8) x Acu ( pP 



(1 + cj) 2 VPS 



substituting this into (IHTj) . we have 



"cq 



1 ± 



4cj /p F 



1/2' 
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p- 



so that 



eq 



P I P 



1/2 



-3/4 



1 ± 



32cj 



(1 + w) 



The radiation energy density p|? today is given by 



30 



-7T 



! iV(Tp)Tp\ 



(85) 



2 1 



1/2 ' 



-3/4 



(86) 



(87) 



where as shown in the Appendix A, N (Tp) = || in the present model (for neutrino filled 
universe). The observational value of T P = 2.7° K yields 



pP = 3.8 x 10 
The total energy density today is 



-51 



(GeV)' 



(88) 



p p = aHl « 4 x 10~ 47 (GeV) 4 , (89) 

where we have used Hp = 72 km s" 1 Mpc -1 . This shows that the universe is matter- 
dominated at present with w p p . This feature is not permanent in the model under 
consideration, since a second era of radiation-dominance starts when 



a — a. 



,ii 

cq 



1 p^il + u) 



3/2 



32n> /pP 
(1+cu) 2 \p\ 



1/2 ' 



(90) 



>/8 2 (pP) 3/2 

With the ratio Pr/Pm ~ 1 x 10" 4 one can approximate the expressions for a eq , T eq and 
a p , obtaining 



l/2 (Af) l/2 



V8 



" ^(l+c) 1 / 2 



1/2 x 2 ' 



(91) 



and 



and 



T ~ 1 T 1 + 
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(1+uj) 



1/2 



a 



1/2 



(l+^x^IxlO 44 (GeV)- 1 



From these we deduce that the 1/a 4 term in p T is dominant at a = a eq , since 



(92) 



cuaz 



~eq 



x - x 3 x 10 7 , 
4 



UJ 



(93) 



for all uj. This implies that one has 1/a 4 dominance throughout a± < a < a cq except perhaps 
for the matter generation periods. Dominance of p T by this term extends well beyond a eq to 
values of a such that a ss 10 _1 cj 1 / 2 a p 3> a cq - However, to decide which of the two terms in 
p r , if any, is dominant at present would depend upon the value of the parameter uj. 

One notes that with 1/a 4 dominance of p r one approximately has aT constant as in the 
standard model. However this does not imply that the entropy is constant under these con- 
ditions, since the variation of aT over the whole range amounts to considerable generation 
of entropy. In the present model there does not exist an entropy problem in any case since 
the entropy is initially zero. In the pure radiation era one has 



with 



and 



we obtain 



so that 



dS_ 

da 



a 

r 



p = — 7r 2 jVT 4 = —Aa 2 - al) , 
H 30 2a 2 V 0J ' 



T 



15a \ 1/4 (I \ , 2 2 ,i/4 



dS fn 2 N\ 1/4 



da \ 15a 



aa 



a z — a 



2\l/4' 

o) 



(94) 



(95) 



(96) 



(97) 



c /1\ 3/4 4 /A 3 V /4 , 2 i 
3= {s) U"30-J ( ° ] 



\ .'- 2 .\n ;i \' ; ,. 2 2 ,3/4 



3 V 30 

where for a > 02, i.e. after matter generation, the entropy will be 



(98) 
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s 



3/4 



4 /vr 2 iV« 3 ^ 1/4 



(d 2 + cud 2 ) 



2\3/4 



+ constant. 



3 V 30 

This gives for the total entropy generated during the era a > a 2 the expression 



(99) 



S (dp) - S (a 2 ) : 
which means that 

5 (a p ) - 5 (a 2 ) 



3/4 



4 / 7T 2 iVa 3 ^ ' 1 

3 1 30 



((i + c) 3 / 4 -M 3/4 ) 



3/4 



3/4 



x 2(l+cu) 3/4 ((l + cu) 0/ - ur J i x 10 



x3/4 



\3/4 



03 



(100) 



It is thus seen that a lot of entropy is produced since the end of the matter generation 
era. We remark that although the parameter u occurs in all expressions for a eq , T eq , a p and 
S (dp) , it is possible to eliminate it and obtain relations between these quantities. One may 
also obtain bounds on these quantities that hold for all u, such as for example, 

d p > 5 x 10 43 (GeVr\ (101) 
This bound leads to an upper bound on the present value of the cosmo logical constant, 



A p < 6 x 10~ 8S (GeV) 



(102) 



which is well within the upper limit of 10 (GeV) placed on A p from cosmo logical obser- 



vation 
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We now consider the time variation of a for a > a 2 . The field equation (Tj~5 



-) =«"Mp r +p m ] + a~ 1 A(i)--^. 
d / 2a z 



or 



a- 1 El u fa 



2d 2 



may be written in the form 



da 
~dt 



1/2 



2(l + o;) \d p 
Let (d p /d) = x and dx = d (d p /d) , therefore 



UJ 



2(l + o;) 



(103) 
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under the approximate boundary condition t = t 2 ~ 0, when a = a 2 ~ 0. For the present 
age t p of the matter-dominated universe one obtains, for all lu, 



( a V f 1 X< ^ X ( \ ( a \ f 1 X 1 / 2 dx fir\c\ 

\pj) Jo Wp-)x + l/2]^ ~ tp {U) ~ Jo [(pVp?) + (1/2) x] 1/2 ' ( } 

Another parameter which is almost u- independent is the present value q p of the decel- 
eration parameter defined by q = —ad /a . One obtains for this value 

ad f d\ /a\ 2 . . 

q --7 = -{a)(a) ■ (106) 



Substituting for ( | j and (|) from the above we obtain 



2 2/ VPm. 

so that g p » | in this model. This is nearly the same value as in the standard model with 
k — 0. 

We now consider the period of generation of rest-mass and in particular show that in this 
model the pressure must have been negative during part of this period. 

C. Period of matter generation 

Consider Eq. ([16]) with A given by (J3 



dE + pda 3 = ada, 

where we have written E for the total energy pa 3 . This equation is valid for all a. Integrating 
it between a = ao and a = a 2 , the end of the matter generation period, one obtains: 
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3 / pa 2 da = a (a 2 — ao) — E 2 



(10£ 



ao 



where E 2 = E (02) and we used the initial condition E (ao) = 0. Now, 



E 2 = E m (a p ) + E T (a 2 ) 



(109) 



since we assume that the rest-mass contribution E m to the total energy remains constant 
for all a>a 2 . Using Eqs. (|6l]h (|57|l and (|55|l, Eq. ffT08l) becomes 



= -En 



+ (a 2 ) 



Pm« P + Pr (fl 2 ) «2 



aa2 



^1(^1+1 

a 2 



2 a 2 ' 



(110) 



Substituting Eq. ffTTUD into ffT08l) we obtain 



3 / pa da = a (a 2 — a ) — E 2 



ao 



(1 + w) 



r m 

p" , , 

/Or / V fl 2 



— — - -aa . (Ill) 
2 a 2 



The right-hand side of fillip is clearly negative since > p P , thus the integral on the 
left-hand side of this equation must be negative. But 



pa 2 da>0 (112) 



indicating that the pressure must have been negative during part of the region of matter 
generation. 

As mentioned before, the region of matter generation is a phase transition period between 
region of known and different equations of state. It now appears that the equation of 
state associated with the creation of rest-mass is characterized by negative pressure. Ideas 
on phase transitions in the early universe from unified gauge field theories indicate that 
phase transitions are expected to have occurred at least twice; at T ~ 10 15 GeV (GUT 
phase transition) and at T ~ 10 2 GeV (electro- weak phase transition). If the maximum 
temperature is of order 10 19 GeV, this gives a\ < 10~ 10 (GeV) 1 so that t\ < 10~ 34 s. 

We now observe that the period of matter generation separates the pure radiation regime 
in which a > and matter-and-radiation regime in which a < as may be seen as 
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This can be written as 



which means that 



and that, for a > a<i 



(X All , . Acff'cctivC 

— = — ^effective + Sp) H . 

a 3 3 



a 
a 



a 



(c + 3p) 



1 



ad = 1 — — (p + 3p) a 2 , 
2a 



— (p + 3p)a 2 



i + -% + ^a 2 >l. 



2 a 



(113) 



(114) 



(115) 



It thus appears that the creation of rest-mass results in the reversal of the sign of a. 
Thus one might say that the presence of rest-mass causes the deceleration of the expansion 
of the universe. The rate of generation of energy is given by the equation 



dE 
da 



a — 3pa^ 



(116) 



This shows that the maximum rate of energy generation occurs within the negative pressure 
interval after the pressure has attained its maximum negative value. 



V. DISCUSSION 



Motivated by previous studies of the back-reaction effect of quantum fields in the Einstein 
static model for the universe we have tried to construct a working model for a dynamic 
universe. The background geometry is assumed to be that of the closed Robertson- Walker 
model. We assumed the existence of a non-zero cosmological constant for start and have 
solved the Einstein field equations accordingly. Quantum field are shown to produce a non- 
zero vacuum energy density which will be a source for a critical density dynamic universe 
that will continue expanding without limit. The non-zero value of the vacuum expectation 
value of the energy-momentum tensor of the quantum field causes the universe to have a non- 
singular start, though the radius is of the size of Planck length. The critical density universe 
was originally considered in Friedmann cosmology to indicate a spatially flat universe, but 
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since we have introduced a non-zero cosmological constant here then it is legitimate not to 
consider the universe to be exactly flat but to be nearly flat. This might explain the recent 
analysis of the CMB measurements which showed that the universe is nearly flat rather 
than being exactly flat. The subsequent development of the universe is controlled by the 
back-reaction effects. Our model is free from all the standard big bang model problems 
and is non-singular. However, it is found that there should be a continuos creation of 
matter and energy in this universe at a rate lower than that proposed by the steady-state 
theory. The rate of matter/energy creation is proportional to the radius of the universe and 
consequently the overall density of universe drops as 1/a 2 . This creation of matter/energy 
can be explained by the fact of the conversion of the energy contained in the cosmological 
constant into particles and radiation. We did not attempt to suggest a mechanism for the 
creation of such energy but we have determined the necessary relations to the cosmological 
constant. 

Figure Caption: 

Fig. 1: The temperature-radius relationship deduced from the back-reaction effect of the 



vacuum energy ploted for different matter fields: the conformally coupled scalar field ( 



), 

the neutrino field (2), the photon field (3) and the minimally coupled field (4) (see ref. [10|). 

Fig. 2: Depicts the contributions of the conformally coupled scalar field (1), the photon 
field (2), the neutrino field (3) and the minimally coupled scalar field (4) to the cosmological 
constant in an Einstein universe at finite temperatures (see ref. [lo|). 

Fig. 3: The temperature-radius relationship according to this model. Note that the 
x-axis does not start from zero but from 0.34. 

Fig. 4: The temperature dependence of the cosmological constant calculated in accor- 
dance with the present model. 
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